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Abstract 

The scaling functions f{ip') and F{y) from the ip'- and y-scahng analyses of inclusive electron 
scattering from nuclei are explored within the coherent density fluctuation model (CDFM). In 
addition to the CDFM formulation in which the local density distribution is used, we introduce a 
new equivalent formulation of the CDFM based on the one-body nucleon momentum distribution 
(NMD). Special attention is paid to the different ways in which the excitation energy of the residual 
system is taken into account in y- and ?/^'-scaling. Both functions, /(V"') and F(y), are calculated 
using different NMD's and are compared with the experimental data for a wide range of nuclei. 
The good description of the data for y < and ip' < (including ip' < —1) makes it possible 
to show the sensitivity of the calculated scaling functions to the peculiarities of the NMD's in 
different regions of momenta. It is concluded that the existing data on the ip'- and y-scaling are 
informative for the NMD's at momenta not larger than 2.0 -i- 2.5 fm^^. The CDFM allows us to 
study simultaneously on the same footing the role of both basic quantities, the momentum and 
density distributions, for the description of scaling and superscaling phenomena in nuclei. 

PACS numbers: 25.30.Fj, 21.60.-n, 21.10.Ft, 24.10.Jv, 21.65. +f 



1 



I. INTRODUCTION 



The inclusive scattering of electrons as weakly interacting probes from the constituents of 
a composite nuclear system is a strong tool in gaining information about the nuclear struc- 
ture. This concerns particularly the studies of such basic quantities of the nuclear ground 
state as the local density and the momentum distributions of the nucleons. As known 111 M 

nnn 

see also UMm) the mean-field approximation (MFA) is unable to describe simultaneously 
these two important nuclear characteristics. This imposes a consistent analysis of the role of 
the nucleon-nucleon correlations using theoretical methods beyond the MFA in the descrip- 
tion of the results of the relevant experiments. It was realized that the nucleon momentum 
distribution (NMD), n{k), which is related to both diagonal and non-diagonal elements of 
the one-body density matrix is much more sensitive to the nucleon correlation effects than 
the density distribution p(r) which is given by its diagonal elements. Thus it is important 
to study these two basic characteristics simultaneously and consistently within the frame- 
work of a given theoretical correlation method analyzing the existing ernpirical data. Such 
a possibility appears in the coherent density fluctuation model (CDFM) P, O, IT 
is related to the delta-function limit of the generator coordinate method (see also 
main aim of the present work is to apply the CDFM to the description of the experimental 
data on the inclusive electron scattering from nuclei, which showed scaling and superscaling 
behavior of properly defined scaling functions, and to gain more information on the NMD 
and the density distributions in nuclei. 

In the beginning of Section IH] we will review briefly the scaling of both first and second 
kind. Scaling of the first kind means that in the asymptotic regime of large transfer momenta 
g = |q| and energy u a properly defined function of both of them F{q, uj) (which is generally 
the ratio between the inclusive cross section and the single-nucleon electromagnetic cross 
section) becomes a function only of a single variable, e.g. y = y{q,uj). This is what is 



I, 8 1 which 
. The 



called y-scaling (see,e.g. 



:tion only oi a single variable, e 

Qfl Q Q, Q H H 0, Q 



). Indeed, for the region y < and 



q > 500 MeV/c this scaling is quite well obeyed. It has been found that the scaling function is 
related to the NMD, and thus some information (though model-dependent) can be obtained 
from the y-scaling analysis. Another scaling variable ip' (related tow) and the corresponding 
'?/''-scaling function f{ip') were defined and considered (see, e.g. 
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20 



21 



within the 



framework of the relativistic Fermi gas (RFG) model. It was found from the studies of the 
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inclusive scattering cross section data that /(V'O shows for ip' < both scahng of the first 
kind (independence of q) and scahng of the second kind (independence of the mass number 
A for a wide range of nuclei from ^He to ^^^Au). This is the so called superscaling jioj ]. 
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2J| . Here we 



The extension of the ■^/''-scaling studies using the RFG model was given in 
would like to emphasize that, as pointed out in 2l|, the actual nuclear dynamical content 
of the superscaling is more complex than that provided by the RFG model. For instance, 
the superscaling behavior of the experimental data for f{ip') has been observed for large 
negative values of ip' (up to ip' ~ —2), while in the RFG model f{ip') = for ip' < —1. This 
imposed the consideration of superscaling in theoretical approaches which go beyond the 
RFG model, i.e. for realistic finite nuclear systems. Such a work was performed using the 
CDFM in j9|. The calculations in the model showed a good quantitative description of the 
superscaling in finite nuclei for negative values of ip', including those smaller than —1. We 
would like to note that the main ingredient of the CDFM (the weight function) was expressed 
and calculated in on the basis of experimentally known charge density distributions p(r) 
'or the ^He, ^^C, ^^Al, ^^Fe and ^^"^Au nuclei. At the same time, however, we started in 
^ the discussion about the relation of /{'ip') with the NMD, n{k), showing implicitly how 
f{ip') can be calculated on the basis of n{k). In Ref. l9] we indicated an alternative path 
for defining the weight function of the CDFM which is built up from a phenomenological or 
a theoretical momentum distribution. In the present paper we give (in Section |n} and use 
(in Section ITTT]l the explicit relationship of flip') with n{k) using the basic scheme of the 
CDFM and showing also how information about n{k) can be extracted from the ip'-scaling 
function. We point out in our theoretical scheme and in our calculations the equivalence of 
the cases when f{ip') is expressed through both density p{r) and momentum distribution 
n{k). In this way both basic quantities are used and can be analyzed simultaneously in the 
studies of the scaling phenomenon. 

Additionally to [9], we present calculations of f{ip') for q = 1560 MeV/c and the com- 
parison with the experimental data from 2^ . 

In the present work we also define the ?/-scaling function F{y) in the CDFM (Section ITT|l 
and present the comparison with the experimental data (taken from jisl . 14|) of our calcu- 



lations of F(y) based on three different NMD's: from the CDFM, from the ^/-scaling (YS) 



studies in Il3, 



14l | and from the parameter-free theoretical approach based on the light-front 



dynamics (LFD) method (Section lllljl . We discuss the sensitivity of the calculated 
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function F{y) to the peculiarities of the different NMD's considered. 

We also estimate the relationship of /{ip') with F{y) and show in Section UTTl the condition 
under which the NMD ncw{f^) extracted from the YS analyses jl^, ll^ can describe the 
empirical data on f{ip')- 

The consideration of the points mentioned above made it possible to estimate approxi- 
mately the region of momenta in n{k) which is mainly responsible for the description of the 
y- and ^'-scaling and how it is related to the experimentally studied regions of the scaling 
variables y and ip'- 

The conclusions of the present work are summarized in Section IIVI 



II. THE THEORETICAL SCHEME 



We start this Section with a brief review of the y- and ^/''-scaling analyses in Subsec- 
tions and respectively. An important point that we emphasize is the way in which 
the excitation of the residual system is taken into account, which is different in each case. 
We discuss the peculiarities of both approaches that are necessary to take into account for 
the development performed in our work within the CDFM f Subsections UTTl and ITTD|1 . 



A. Brief review of the y-scaling 

In this Subsection we will outline the main relationships which concern the ly-scalin g in 
the inclusive electron scattering of high-energy electrons from nuclei (e.g. 

mm 

Q, Q, At large transfer momentum (g > 500 MeV/c) and transfer energy u, 

the scaling function F{q,uj), which is the cross section of the inclusive process divided by 
the elementary probe- constituent cross section, turns out to be a function of only a single 
variable y = y{q,uj). This is the scaling of the first kind. The smallest value of the missing 
momentum p = |p| = Ipat — q| (pat being the momentum of the outgoing nucleon) at the 
smallest value of the missing energy is defined to be y {—y) for u larger (smaller) than its 
value at the quasielastic peak 

Ljc^iq^ + miy/^ -rriN, (1) 
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mjv being the nucleon mass. The condition for the smallest missing energy means that the 
value of the quantity 

g(p) = ^(Ma-i^+p' - ^iMUy+? (2) 

(where Ma_i is generally the excited recoiling systems's mass and M^_]^ is the mass of the 
system in its ground state) must be 

S{p) = 0. (3) 

The quantity S{p) Q characterizes the degree of excitation of the residual system and 
essentially it is the missing energy {Em) minus the separation energy (Es). So, at the 
condition (jHj) E^ = E^ . 

As shown (e.g.|lll[liQ,Q), for q > 500 MeV/c 

F{q,y)'^F{y) = fiy)-Biy), (4) 

where 

/•OO 

f{y) = 2n n{k)kdk (5) 
and n{k) is the conventional NMD function normalized to unity 

1. (6) 



j dkn(k) 



The information on F{y) and, correspondingly, on f{y) can be used to obtain n{k) by: 

1 df{y) 



n{k) 



2ny dy 



(7) 

y\=k 



In Eq. (0]) B{y) is the binding correction which is related to the part of the spectral 
function generated by ground-state correlations and the excitations of the residual system 
(when Ma-1 > ^^d, correspondingly, S{p) > ). 

The problem of the correct account for the binding correction is a longstanding one. 
Only when the excitation energy of the residual system is equal to zero (as in the case of 
the deuteron) B = and then F{y) = f{y). Generally, however, the final system oi A — 1 
nucleons can be left in all possible excited states. Then B{y) and F{y) ^ f{y)- 

In a new y-scaling variable ( yew) "was introduced on the basis of a realistic nuclear 



spectral function as provided by few- and many-body calculations 
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27l |. The use of yew 



leads to B{ycw) = and, consequently, to F{ycw) = fivcw)- The latter is important 



because in this case it becomes possible to obtain information on the NMD directly (using 
Eq- ©) without introducing theoretical binding correction B{y). In this consideration the 
removal energy (whose effects are a source of scaling violation, the other source being the 
final-state interactions) is taken into account in the definition of the scaling variable. So, 
the binding corrections are incorporated into the definition of yew- 

The analysis of empirical data on inclusive electron scattering from nuclei (with A < 56) 
that the following form of f{y) gives a very good agreement with the data: 



showed 



HE 



fiy) 



Ci exp{-a^y^) 



+ C2expi-b\y\){l + 



by 



(8) 



+ Y 1 + y^/a"^ J 

where the first term describes the small y-behavior and the second term dominates large y. 
From Eq. ((7j) one can obtain the following form of n{k): 



n{k) = UMFAik) + Ucorrik), 

where the mean-field part nMFA{k) of the NMD (for k ^ 2 fm~^) is: 



nMFA{k) 



TC 



1 + a^(a^ + r 



exp(— a^/E^^ 
(a2 + A;2)2 



(9) 



(10) 



while the high-momentum components of n{k) which contain nucleon correlation effects are 
given by: 



rir 



ik) 



exp{—bk) 
27r(l + P/a2) 



3 + bk + 



(11) 



30ut the n{k) from the y-scaling analysis 



Further in our work we will use the information a 
[Eqs. fl^- (fTT|) ]. The values of the parameters |l^, e.g. in the case of interest for the 
s^Fe nucleus are: b = 1.1838 fm, Ci = 0.30 fm'^ C2 = 0.11838 fm, a = 0.710 fm-^ and 
a = 0.908 fm. 



B. The V' -scaling variable and the V' -scahng function in the relativistic Fermi gas 
model. Relation between the y- and ^'-scaling variables 

In this Subsection we will review briefly the scaling in the framework of the RFG model 
19,0, 21, Q- 

This will be necessary for our consideration of the scaling in the present 
work within the CDFM in the next Subsections III CI and III Dl The ?/-scaling variable in the 
RFG has the form 

Vrfg = niN I AW 1 + - - k] , (12) 



where 

K = q/2mN, A = uj/2mN, r = |(5^|/4m^ = — 

are the dimensionless versions of q, uj and the squared four- momentum \Q'^\ 
22| | a new scahng variable if) was introduced by 

1 A-r 



(13) 



In 



2G 



21 



where 



(14) 



(15) 



= y 1 + rjp — 1 and rjp = kp/m 

are the dimensionless Fermi kinetic energy and Fermi momentum, respectively. 

In order to include, at least partially, the missing e nerg y dependence in the scaling vari- 
able, a shift of the energy cu is introduced in the RFG j2l|: 



LU' = UJ-E 



shifti 



(16) 



where Eghift is chosen empirically (in practice it is from 15 to 25 MeV) and thus can take 
values other than the separation energy Eg. The corresponding A and r become 



A' = Co''/2mAr, T 



A 



l2 



(17) 



This procedure aims to account for the effects of both binding in the initial state and 



2l| that the corresponding new version 



interaction strength in the final state. It is shown in 
of the ■j/'-scaling variable {ip') has the following relation to the y-scaling variable: 

Z/oo(A = A') 



^' = V^[A - A'] 



A 



kp 

UJ 



11 yoo(A = AO 
^ + r + 4;^2^" kp 

ijj - E^ 



(18) 



2mj\f 2m jv 

In (11 8j) Uoc, is the y-scaling variable in the limit where M^_^ — »• oo. kp is the Fermi 
momentum which is a free parameter in the RFG model, taking values from 1.115 fm~^ for 
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C to 1.216 fm-i for ^^'^An 



21[. As shown in 2l|, Eq. ()18|) contains an important average 



dependence on the quantity £{p) (0) (i.e. on the missing energy) which is reflected in the 
quadratic dependence of ip' on the ?/-scaling variable. 



Finally, in 



21 



a dimensionless scaling function is introduced within the RFG model 



fRFOm = kpFRpG^. 



(19) 



The careful analysis of the experimental data on inclusive electron scattering [2l|, |22] shows 
that the RFG model contains scaling of the first kind (/ or F are not dependent on q 
at high-momentum transfer and depend only on ip') but also that f{ip') is independent of 
kp to leading order in rip, thus showing no dependence on the mass number A (scaling of 
the second kind). In the RFG both kinds of scaling occur and this phenomenon is called 
superscaling. 



To finalize this Subsection we give the analytical form of fnFG obtained in 
which will be used in this work: 



IS 



2G 



21 



2x1 



Vf 



(20) 



4- ■ ■ Vf 

Here we note that due to the 0-function in (j2(jp . the function f{ip') is equal to zero at 
ip' < —1 and ip' > 1. As can be seen in Fig. 1 of Ref. |9||, this is not in accordance with the 
experimental data and justifies the attempt in jsl, as well as the development made in the 
present work, to consider the superscaling in realistic systems beyond the RFG model. 



C. Theoretical scheme of the CDFM and the ip' -scaling function in the model 

The CDFM was suggested and developed in 

mm 

. It was deduced from the delta- 
function limit of the generator coordinate method |28|. The model was applied to the study 
of the superscaling phenomenon in In this Subsection we continue the development of 
the model aiming its applications to the studies of the NMD from the analyses of the y- and 



T/^'-scaling in inclusive electron scattering from nuclei. We will start with t 
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le expressions 
). They are 



of the Wigner distribution function (WDF) in the CDFM W{r, p) (e.g. 
based on two representations of the WDF for a piece of nuclear matter which contains all 
A nucleons distributed homogeneously in a sphere with radius R, with density 

PoW = ^ (21) 

and Fermi momentum 

kF = kF{R) = (^Po(i?)) - ^, « = (^) ^ 1.52AV3. (22) 
The first form of the WDF is: 

Wniv, p) = j^,m - l^rnkFiR) - IpI). (23) 



The second form of the WDF for such a piece of nuclear matter can be written as 

W^I,(r,p) = (^©(^^ - IPl)0(^ - IpD- (24) 

In the CDFM the WDF, as well as the corresponding one-body density matrix (ODM) can 
be written as superpositions of WDF's (ODM's) from Eqs. (j^^ and (j^^ in coordinate and 
momentum space, respectively: 

iy(r,p)=y^ dR\F{R)\^WR{v,^) = —--j^ dR\F{R)\^Q{R-\v\)Q{kF{R)-\v\) (25) 



and 



/■oo A roo Q, 

W{v,p)= dkF\G{kF)\'WT:^{r,p) = -—J dkF\G{kF)\'Q{kp-\p\)ei^-\r\). 
Jo ^ [Ztt)-^ jo kp 

(26) 

The relationship between both |Fp and |Gp functions is: 

|G'(MP = A|i^(f )r- (27) 
kp i^F 

Using the basic relationships of the density and momentum distributions with the WDF: 

p(r)= JdpWiv,p), (28) 

n(p) = J drW{r,p), (29) 
one can obtain the corresponding expressions for p(r) and n{p) using the WDF from Eq. (j23): 

/•CO Q A 

p(r)=y^ dR\F{R)\'^Q{R-\r\), (30) 
nip) = :^ r^'dR\FiR)\'R', (31) 



Svr^ Jo 

and, equivalently, using the WDF from Eq. (pUj): 



p(r) = — / dkp\G{kp)\^Q{- - kp)kp\ (32) 
Jvr^ JO r 

n(p) = / rfA;^|G(MI'^30(^^ " IpD' (33) 
JO 47rA^r, 



Ank 

both normalized to the mass number 



J p(r)dr = A, J n{'k)d'k = A (34) 



when both weight functions are normahzed to unity: 



dR\F{R)f 



(35) 



One can see from Eqs. (jHUj) . (jHH) and (jHHj) the symmetry of the expressions for p(r) 

and n{p) as integrals in the coordinate and momentum space. 

A convenient approach to obtain the weight functions F{R) and G{kF) is to use a known 
(experimental or theoretical) density distribution p(r) and/or the momentum distribution 
n{p) for a given nucleus. For |F(i?)p one can obtain from Eqs. and (jHH): 

1 dp{r) 



(at dp/dr < 0) and 



\FiR)r 



Po{R) dr 
Stt^ a dn{p) 



r=R 



p=a/R 



2 i?5 dp 
(at dn/dp < 0). 

The expressions for |G(A;j7)p can be obtained from Eqs. (jH^ and 

Stt^ a (ip(r) 



(36) 



(37) 



(at dp/dr < 0) and 



(at dn/dp < 0) with 



|G(fci 



|G(fci 



2 kp' 



1 



dn{p) 



r=a/kf 



no{k. 



rioikp) dp 

, 3A 
V — — 



p=k[ 



(38) 



(39) 



(40) 



In order to introduce the scaling function within the CDFM, we assume that the scaling 
function for a finite nucleus fXip') can be defined and obtained by means of the weight 
function |F(i?) p (and \G{kF)\'^) weighting the scaling function for the RFC model depending 
on the scaling variable ip'^^ ifRFGi"^' = V^r); Eq. (j^ ). corresponding to a given density 
Po{R) (PT|) and Fermi momentum kpiR) ^ (or corresponding to a given density in the 
momentum space no{kp) (HO}). 

One can write the scaling variable ipR the form j^: 

p{y) p{y)R 



^Riy) 



kp{R) 



a 



(41) 
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where 



with 



piv) 



y(l + c?/), y>0 
\y\{l-c\y\), y < 0, < l/2c 



1 



1 



2mN V 4/?^ 
Also a more convenient notation can be used: 



kp p{y) _ kp 



(42) 



(43) 



(44) 



kF{R) kp kpiR) 

Using the 9-function in Eq. the scahng function for a finite nucleus can be defined by 
the following expressions: 

ra/{kp\-4,'\) 



dR\F{R)\^fRPG{R.^') 



(45) 



with 



fuFciR,^') 



'kFRW\ 



a 



RmN\^ fkpRlfl 



a 



a 



X 



2 + 



/ a 



\Rm 



N 



2il + 



/ a 



\Rm 



N 



(46) 



and also, equivalently, by 



dkp\GikF)\'fnFGikF,^lj') 



(47) 



with 



fRFGikp,^') 



'kpWl 



ki 



kf 



- 2 



1 + 



(48) 



In this way in the CDFM the scaling function f{ip') is an infinite superposition of the RFG 
scaling functions fRFciR^i'') (or fRFcikp^ilj'))- 

In Eqs. ()45 |) -()48 |) the momentum /cj? is not a free fitting parameter for different nuclei, as 
in the RFG model, but can be calculated consistently in the CDFM for each nucleus (see 
fj36 |l - (|39|l ) using the expression 



kF= dRkF{R)\F{R)\^ = a dR-\F{R)\' 
Jo Jo R 



47r(97r)i/3 
3^2/3 



/•OO 

/ dRp(R)R (49) 

JO 
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when the condition 



is fulfilled and, equivalently, 



when the condition 



3A Jo 



1677 /■°° — _ _ 3 

dkFn{kF)kF , 







(50) 



(51) 



_lim 

kp—*oo 



(52) 

and (jSH) can be used 



n{kF)kF 

is fulfilled. Generally, Eqs. (jSUj) and are fulfilled, so the Eqs. 
to calculate fc^ in most of the cases of interest. 

The integration in and (jTfjl . using Eqs. leads to the following expressions 

for fiij'): 



and 



An r»/(kF\-<P'\) 



A Jo 



Air r°° 



dRp{R) 



R'fHFG{R,^P') + 



B^dfRFciRAn 
3 dR 



A JkFW\ 



dkFn{kF) 



kF fRFcikFj-ip') + 



kF^ dfRFG{kF,i^') 



(53) 



(54) 



0, 



(55) 



3 dkF 

the latter at 

_lim n{kF)kF'^ 

fci?— >oo 

where the functions fuFciRji'') and fRFcikFj'ip') are given by Eqs. (jlFIjl and (jlKjl . respec- 
tively. We emphasize the symmetry in both Eqs. and ()54|) . We also note that the 
CDFM scaling function f{ip') is symmetric at the change of ip' to 

The scaling function f{ip') can be calculated using Eqs. and by means of: i) its 
relationship to the density distribution p(r), and ii) from the relationship to the NMD n{p). 
Both quantities (p and n) can be taken from empirical data or from theoretical calculations. 
In the CDFM they are consistently related because they are based on the WDF of the 
model (Eqs. ()25|) and (j2SI))- Using experimentally known density distributions p(r) for a 
given nucleus one can calculate the weight functions |Fp (Eq. (jHSI)) or \G\'^ (Eq. (jHHj) ) and 
by means of them to calculate n{p) in the CDFM (by Eqs. (jHH) or (jH^ . respectively). 

From Eq. one can estimate the possibility to obtain information about the NMD 
from the empirical data on the scaling function f{ip'). If we keep only the main term of the 
function fRFcikF^i^') 



fRFG{kF,i>') 



kF^ 



(56) 
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and its derivative 



then: 



In Eq. ^ 



dkf 



2 k 



3 ' 



f — — — 2 

f{^lj') ~ Stt / dkFn{kF)kp 



3 If' 



n('kF)d\<.F = 1. 



(57) 



(58) 



(59) 



'L (k I I )^ 

Neglecting the second term in the bracket in ()58|1 (because — ^ 1) one obtains: 

3 kp 



f°° — — — 2 
f{ip') ~ Stt / dkFn{kF)kF ■ 



^'1 

Taking the derivative on from both sides of Eq. ()60|) leads to: 

1 df{i,') 



(60) 



n{p) 



(61) 



\i>'\=p/kf 



37rp^kF d{\ip'\) 

Eq. (|F)T|l can give approximately information on the NMD n{p). If one keeps the second 
term in the bracket under the integral in (j^Hj) . then more complicate equation results: 



dikpm 



= —2'Kp n{p) — 27rp j dk n{k ). 

kp\'il>'\=p ^ 



(62) 



D. y-scaling function in the CDFM and the relationship between the y- and 
scaling functions in the model 

In order to define the y-scaling function F{y) (with dimensions) in the CDFM and to 
establish the relationship between the latter and the dimensionless scaling function f{'4)') 



(Eqs. (E 
model \21 



and (|54j) ) we start with the expression which relates both functions in the RFG 

fB.FG(kF,^'iy)) 



RFG 



(y) 



kf 



(63) 



In analogy with the definition of fiip') in the CDFM, we introduce the function F{y) in a 
finite system as a superposition of RFG y-scaling functions FRpciy) (j^ : 



F{y) 



kFW\=\p{y)\ 



dkp\G{kF)\^ 



fRFG{kFW\ = \p{y)lkF) 

kp 



(64) 



13 



where the function Jrfg has the form ()48p and kplip'l 
Using Eq. ^ for |G'(¥f)P, we obtain from Eq. 



Fiy) 



Air r°° 



with p{y) given by Eq. 
2dfRFGi\piy)\,kF) ] 



(65) 



3 Jkp\tp'\ = \p{y)\ ^ ^ i ^' ' ^ QJ^p J 

Keeping in Eq. the main terms of the function /rfg [Eq. (fHHj) ] and of its derivative 
dfRFc/dkp [Eq. (j37j) ]. one obtains the scahng function F{y) in the form: 



F{y) ~ 27r / dkkn{k). 
In Eqs. fl65|) and ()66|) the normahzation of n{k) is: 

y n{k)dk = 1. 
For the cases of interest when y < {\y\ < l/(2c)): 

/■oo 
c/A; A;n(fc). 
,y|(i-c|j/|) 

If the T/^'-scaling variable is not a quadratic function of y as in |2ll . l2z 
it is a linear one (i.e. kpli^'l = \y\), then 



(66) 



(67) 



(68) 

(see Eq. ((THD) but 



POO 

F{y) ^2tx dkkn{k). (69) 

Eq. flU^ gives the known y-scaling function jl^ ll^ and its relationship with 

the NMD n{k). It can be seen from ()68|) that the use of a more complicated (quadratic) 



dependence of ip' on y 



22| leads to a more complicated ?/-scaling function in the CDFM 



and its relationship with the NMD. 

To finalize this section we elaborate somewhat more on Eq. (j6Up , considering the relation- 
ship between the scaling variables y and ip' and between the corresponding scaling functions. 
Eq. ()60|) for the T/^'-scaling function f{ip') can be rewritten in the following approximate form: 



2 3 

dkpn{kp)kp ~ -kav^n 



kp 



^'1 



If we admit as in the case of the |/-scaling 



10 



11 



dkpn{kp)kp. 



kpWl = \y\ 



then 



(70) 



(71) 



(72) 
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where F{y) is the ?/-scahng function (jU^ and kav can be estimated as in 



K.^{ly\ where (I) = [ dl.^. (73) 



In the case of the ijj' -scaling variable (for y < 0) 

kpm = \y\il-c\y\) (74) 

and we can replace approximately the lower limit of the integration by |?/| which is, however, 
the solution of Eq. 

f{iP') ~ Stt r dkp n(kF)kF^, where \y\ = ^ (l - Jl - Ackpl^'l) . (75) 

J\y\ 2C \ J 

III. RESULTS OF CALCULATIONS AND DISCUSSION 

We begin this Section with calculations of the scaling function /(■?/'') using Eqs. and 
fl54|) for different nuclei within the CDFM for the transfer momentum q = 1560 MeV/c. 
The results for ^He, ^^C, ^^Al and ^^''Au are presented in Fig. ^ and are compared with 



the experimental data from 



values of kp from 
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22| and with the predictions of the RFG model [Eq. ()20|)] with 
22]. Our calculations are performed in addition to those for q = 1000 
and 1650 MeV/c in 2| which were compared with the data from 21]. 

Following the main aim of our work, namely to extract reliable information on the NMD 
from the scaling function, we will consider in detail the consecutive steps of the calculations 
of fiip') in connection with n{k). We note firstly that /{ip') is calculated in the CDFM using 
Eq. where the density distribution p(r) is taken to be of a Fermi-type with parameter 
values obtained from the experimental data on elastic electron scattering from nuclei and 
muonic atoms. For ^He and ^^C we used a symmetrized Fermi-type density distribution [2^ 
with half-radius (-R1/2) and diffuseness (6) parameters: R1/2 = 1.710 fm, b = 0.290 fm for 
^He and R1/2 = 2.470 fm, b = 0.420 fm for ^^C. These values of the parameters lead to charge 
rms radii equal to 1.710 fm for ^He and 2.47 fm for ^^C which coincide with the experimental 



ones 
from 



30|. For the ^^Al nucleus the values of R1/2 = 3.070 fm and b = 0.519 fm ar_e taken 



30[. For ^^^Au the parameter values are R1/2 = 6.419 fm |31| and b = 1.0 fm |9|. The 



H and 6= 1.0 fm Q. 



necessity to use the latter value of b for ^^^Au instead of 6 = 0.449 fm 



3 



was discussed in 



P]. This is an ad- hoc procedure in order to obtain high- momentum components of the NMD 
n{k) (using and (jSSl)) which are similar to those in light and medium nuclei. This was 
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FIG. 1: Scaling function /{^p') in the CDFM (solid line) at q = 1560 MeV/c for ^He, ^^c, 27 ai, 
and ^^''Au. The results are obtained using Eqs. and equivalently by Eqs. with 

n(p) from the CDFM [Eq. The experimental data from 22] are given by the shaded area. 

The RFG result [Eq. 1)20(1 ] is shown by dotted line. 

necessary to be done due to the particular A-dependence of n{k) in the CDFM resulting in 
lower tails of n{k) at k > 2 fm~^ for the heaviest nuclei which has to be improved. Also for 
^^Fe a better agreement with ib' scaling data at q = 1000 MeV/c is obtained for b = 0.7 fm 

n 

instead of 6 = 0.558 fm |3(| and this result will be shown later on. 

As can be seen, the CDFM results for the scaling function f{tp') agree well with the 
experimental data taken from inclusive electron scattering 0|. This is so even in the 
interval ip' < —1 for all nuclei considered, in contrast to the results of the RFG model where 
fRFci'^') = for tp' < —1. Here we emphasize that our scaling function f{ip') is obtained 
using the experimental information on the density distribution. At the same time, however, 
fiip') is related to the NMD n(p), as can be seen from Eq. (j54|l . We note that Eqs. (j53| 
and ()54|) are equivalent when we calculate in the CDFM the NMD n{p) consistently using 
Eq. ()31|) . where the weight function is calculated using the derivative of the density 
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distribution p(r) (Eq. ()36|) ). The same consistency exists in the calculations of the CDFM 
Fermi momentum kp (which is used in the calculations of /{ip') from Eqs. (j^Hj) and 
It is calculated by means of Eq. and, equivalently, by Eq. ^T\i where the CDFM result 
for n{p) is used. The calculated values of kp in the CDFM are: 1.201 fm^^ for ^He, 1.200 
fm-i for 12c, 1.267 fm'^ for ^^Al, 1.270 fm'^ and 1.200 fm-^ for ^^^An. 

In Fig. 121 we give the results of the calculations of the NMD n{k) within the CDFM 
for ^He, i^C, ^^Al, ^^Fe and ^^''Au using Eqs. (jHT|) . (jHUI) and Fermi-density distribution p 
with parameter values mentioned above (with b = 1.0 fm for ^^''Au). The normalization 
is /n(k)rfk = 1. The nik) from CDFM for the nuclei considered are with similar tails at 
k ^ 1.5 fm^i, so they are combined and presented by a shaded area. As can be expected, this 
similarity of the high-momentum components of n{k) leads to the superscaling phenomenon. 
In our work there is an explicit relation of the scaling function to the NMD in finite nuclear 
systems [Eq. (jHH)]- As can be seen, when the latter is calculated in a realistic nuclear 
model accounting for nucleon correlations beyond the MFA, a reasonable explanation of the 
superscaling behavior of the scaling function for ^z^' < — 1 is achieved. 

In Fig. 121 we give also: i) the "^/-sc aling data" for n{k) in ^He, ^^C and ^^Fe obtained 
from analyses of (e, e') cross sections in 12| on the basis of the y-scaling theoretical scheme, 
ii) n{k) calculated within the MFA using Woods-Saxon single-particle wave functions for 
^^Fe^Jii) the NMD for ^^Fe [Eq. ©] extracted from the more recent ^/-scaling analyses in 
0,^^. We give also the corresponding contributions to n{k), namely the mean-field one 
nMFAik) [Eq. (fTn|l ] and the high-momentum c omp onent ricorrik) [Eq. (|TT|) ]. iv) the NMD 
n{k), e.g. for ^^Fe obtained within an approach |25l[ ba sed on the NMD in the deuteron from 



the light-front dynamics (LED) method (e.g. |32l . 33| and references therein). In 25 1 n{k) 
was written within the natural-orbital representation Is^ as a sum of hole-state {n^{k)) and 
particle-state {pP{k)) contributions 

n{k) = Na [n\k) + nP{k)] . (76) 

In dZHl) 

n\k) = ^'2(2j + l)XnijC{k)\Rni,{k)\', (77) 

nlj 

where F.L. denotes the Fermi-level, and 

C{k) = (78) 

{2nf^k^ + m% 
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FIG. 2: Nucleoli momentum distribution n(k) from: i) CDFM: the calculated results using Eqs. (|3H1 
and (15^ for ^He, ^^C, ^'^Al, ^^Fe and ^^'^Au are combined by shaded area {tlcofm)] ii) "y-scaling 



data" m 



analyses 



nven by open squares, circles and triangles for ^He, ^^C and ^^Fe, respectively; iii) 



l3,Q for s^Fe [Eq. ®] {new), umfa [Eq. (CHI)], ncow [Eq. (CU)]; iv) LFD approach |25| 
[Eqs. (Iinil-llZni), (IHH), fo^ ^^Fe {ulfd); v) MFA calculations using Woods-Saxon single-particle 
wave functions for ^^Fe {nws)- 

In Eq. (f77j) Xnij are the natural occupation numbers (which for the hole-states are close to 
unity and were set to be equal to unity in *25] with good approximation) and the hole-state 
natural orbitals Rnij{k) are replaced by single-particle wave functions from the MFA. In 
Woods-Saxon single-particle wave functions were used for protons and neutrons. The use of 
other s.p. wave functions (e.g. from Hartree-Fock-Bogolyubov calculations) leads to similar 
results. 
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The normalization factor has the form: 

-1 

(79) 



Na = Un l^dqq 



F.L. A 

J2 2(2j + l)Xni,Ciq)\R^i,{q)\^ + -n,{q) 



nlj 2 



The well-known facts: i) that the high-momentum components of n{k) caused by short- 
range and tensor correlations are almost completely determined by the contributions of 
the particle-state natural orbitals (e.g. [2^), and ii) the approximate equality of the high- 
momentum tails of n{k)/A for all nuclei which are rescaled version of the NMD in the 
deuteron nd{k) (36j : 

TiAik) ~ aAUdik) (80) 

(where a a is a constant) made it possible to assume in 2^ that n^{k) is related to the 
high-momentum component n^{k) of the deuteron: 

n^k) = ^n,{k). (81) 
In and (|HT|) n^i^k) is expressed by an angle-averaged function 



n,{k)=C{k){l-z^)fi{k). (82) 

In Eq. (jH^ z = cos(k, n), n being a unit vector along the 3- vector (^ component of the 
four- vector u which determines the position of the light-front surface [22l0]- The function 
f^lk) is one of the six scalar functions /i_6(A;^, n-k) which are the components of the deuteron 
total wave function \I'(k, n). The component /s exceeds sufficiently other /-components for 
k > 2 ^ 2.5 fm~^ and is the main contribution to the high-momentum component of nd{k), 
incorporating the main part of the short-range features of the nucleon-nucleon interaction. 

As can be seen in Fig. |2] the calculated LFD momentum distributions are in a good 
agreement with the "y-scaling data" for ^He, ^^C and ^^Fe from jl^, including the high- 
momentum region. We emphasize that n{k) calculated in the LFD method does not contain 
any free parameters. 

The comparison of the NMD's from the CDFM, LFD and from the y-scaling analysis 
(YS) [Eqs. (jm)- (fTT|) ] shows their similarity for momenta k ^ 1.5 fm~^ (i.e. in the region 
where the MFA is a good approximation). It also shows their quite different decreasing 
slopes for k > 1.5 fm~^, where the effects of nucleon correlations dominate. In the rest of 
this Section we will try to consider in more details the questions concerning the reliability of 
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the information about the NMD obtained from the y- and ?/^'-scahng analyses. This concerns 
the sensitivity of such analyses and also the identification of the intervals of momenta in 
which n{k) can be obtained with more reliability from the experimental data and from the 
y- and ^/''-scaling studies. First of all, we emphasize that the approaches considered to 
obtain experimental information on n{k) are strongly model dependent. In this respect we 
note various ways to introduce the scaling variables, e.g. the different y-scaling variables 



l^llMUJl; the different ^/'-scaling variables (e.g. in [l9 



IJ], the different ^/'-scaling variables (e.g. in (l9|,|2l|), as well as the corresponding y- 
and ^/'-scaling functions. Even so, however, it is nevertheless worth considering in more detail 
the model-dependent empirical information about the NMD coming from y- and T/'-scaling 
analyses. We emphasize that our consideration is based on f{ip') and the y-scaling function 
F{y) within the CDFM, as well as on the relationship between both of them discussed in 
Subsection IIIDl 

Firstly, we give the results of the calculations of the ?/-scaling function Fiy) obtained in 
the CDFM [Eq. ^] using different NMD n{k): i) from the CDFM (Eqs. dSU and (EHl)), ii) 
from the yew scaling approach jl^, 14| [Eqs. (Pl- (fTT|l ]. and iii) from the approach which 
uses the results of the LED method [Eqs. (|7n|l - (f7n|) . (jHH), (jH^ ]. In Fig. El they are compared 
with the yew scaling data for Fiy) for the ^He and ^^Fe nuclei taken from As can 

be seen from Fig.Ol there is a general agreement with the data for all the NMD's considered. 
At first thought this can be surprising knowing the different behavior of n{k)^s for larger k 
that are seen in Fig. [2l The reasons for the relative similarity of the results for F{y) using 
different n(A;)'s follows. 

i) The tp'- scaling variable is a quadratic function of y but not a linear one [Eq. (fT^ ]. 
n accordance with this, the lower limit of the integral [Eq. ()65p] for F{y) is not \y\ as in 
13, [1^, but \y\{l — c\y\) (see for a comparison Eqs. and (jSHI)); h) Due to the steep 



slope rates of decreasing of the NMD's for large momenta, the main contribution to the 
integral (j65|) (and to the estimation (j68|) ) comes from momenta which are not much larger 
than the lower limit of the integration. In this way, the very high-momentum components of 
n{k) do not play so important a role (in the integral in (jHSI)), at least for momenta studied 
so far y > —700 MeV/c. We give some numerical estimations: for example, for y = —300 
MeV/c instead of integrating from \y\ = 300 MeV/c = 1.52 fm~^ (as in ()69p). in F{y) in 
()68|) the integration starts from \y\{l — c\y\) = 1.19 fm~^, for y = —600 MeV/c instead of 
integrating from \y\ = 600 MeV/c = 3.04 fm~^ the lower limit of the integral in (jHH|l is 
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FIG. 3: The y-scaling function F{y) for ^He and ^^Fe calculated in the CDFM [Eq. (|65)) ] (solid and 
thick dashed lines), from the j/nw -scalin^ approach Q, Q [Eqs. &-^] for ^^Fe (dash-dotted 
line) and from the approach 25] within the LFD method [Eqs. H76p - (|79p . H81() . (|82l)1 fthin solid 



and dashed lines). The results are compared with the j/cpiz-scaling data taken from 

\y\{l — c\y\) = 1.71 fm~^. This means that the main contribution to F{y) from n{k) is for 
momenta k ^ 2 fm~^. The behavior of F{y) in Fig. IHl reflects that one of n{k). For instance, 
for -400 < y <0 MeV/c the CDFM result for F{y) is higher than those of LFD and YS 
because the values of n{k) from CDFM for k < 1.5 fm~^ are larger than those of n{k) from 
the LFD and the YS. In contrast to this, the values of F{y) for —700 < y < —400 MeV/c in 
the CDFM result are lower than those of the LFD because n{k) from LFD has a higher tail 
than n{k) in the CDFM for k > 1.5 fm~^. Nevertheless, though the tails of n{k) for large k 
are quite different (for k > 1.5 fm^^), the values of F{y) from the different approaches are 
quite close to each other and are in agreement with the existing data. In this way, we can 
conclude from our experience that the existing y-scaling data can give reliable information 
for the NMD for momenta not larger than 1.5-^2.0 fm~^, where the considered n{k) are not 
drastically different from each other. 

One can see from Fig. IHl that the CDFM results for F{y) are in a very good agreement 
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with the data for ^He for y ^ —400 MeV/c, while in the same region the result of the LFD 
agrees very well with the data for ^^Fe. The YS result for F[y) agrees well with the data 
for s^Fe for y > -600 MeV/c. 

It is worth mentioning that in our approach we start from the ■^/''-scaling consideration for 
the function F{y) and this leads to a relatively good description of the y-scaling data on the 
basis of the correct accounting of the relationship between the and y-scaling variables. 
The overall agreement of the theoretical results using the momentum distributions from 
the CDFM, the LFD and the YS with the experimental data for F{y) is related with their 
similarities up to momenta k = 1.5 ^ 2.0 fm~^. 

Our next step is to estimate the ■^/''-scaling function f{ip') [Eqs. (IHU)) and ^7U\\ ] replacing 
the lower limit of the integration kplip'l approximately by |?/| which is, however, a solution of 
(fTIj) . i.e. \y\ ~ 2^ ^ y^l — ickplip'lj, but is not the linear function of \ip'\: \y\ = kplip'l- 
This is done in order to introduce in the relationship of fiip') with F{y) in (fTOj) and ()72|) 
the lower limit in the integral for F{y) to be \y\ (as in the YS) where, however, the correct 
relationship of \y\ with (Eq. (jZH)) is accounted for. In Fig. ID we give the results for 
f{ij') from Eq. JCSI) using the NMD from the YS analysis QQ [Eqs. ©-(HH)] and from 
the approach |25| within the LFD method [Eqs. dHU), ©]• One can see that 

the NMD from the YS analysis [Eqs. ©-(^^j gives a good description of f{ip') for ^^Fe 
in the case of g = 1000 MeV/c for values of ip': — LIO < ip' < (for which y < and 
\y\ < l/(2c) at c = 0.144 fm). The scaling function f{ip') calculated by n{k) from the LFD 
is in agreement with the data for —0.5 ^ 'i/'' < 0, while in the region —1.1 <■?/''< —0.5 
shows a dip in the interval —0.9 < ip' < —0.6. The difference in the behavior of /{ip') in 
these two cases reflects the difference of the momentum distributions of YS and LFD in the 
interval 1.5 ^ k ^ 2.5 fm~^: the n{k) of the LFD has a dip around k ^ 1.7 fm~^ below the 
curve of n{k) from the YS analysis. 



IV. CONCLUSIONS 



The results of the present work can be summarized as follows. 

(i) The main aim of our work was to study the nucleon momentum distributions from the 
experimental data on inclusive electron scattering from nuclei which have shown the phe- 
nomenon of superscaling. For this purpose we made an additional extension of the coherent 
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FIG. 4: The -(/''-scaling function /(^') at ^fFe _and q = 1000 MeV/c calculated from Eq. ((7S|l usin; 
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from: i) the yciy-scaling analysis ^ [Eqs. (|9|)- (|llj) ] (solid line); ii) the approach 
within the LFD method [Eqs. ^] (dashed hne). The CDFM result obtained 

using Eqs. (E I and (jlB)) is given by dotted line. The experimental data given by shaded area are 



taken from 



2l|. 



density fluctuation model in order to express the '?/''-scaling function, fiip'), explicitly in 
terms of the nucleon momentum distribution for realistic finite systems. This development 
is a natural extension of the relativistic Fermi gas model. In this way f{ip') can be expressed 
equivalently by means of both density and momentum distributions. In ^9] our results on 
fiip') were obtained on the basis of the experimental data on the charge densities for a wide 
range of nuclei. In the present work we discuss the properties of n{k) which correspond to 
the results for fiip') obtained in the CDFM. Thus we show how both quantities, the density 
and the momentum distribution, are responsible for the scaling behavior in various nuclei. 

(ii) In addition to the work presented in Ref. [^, the scaling function f{ip') is calculated 
here in the CDFM at g = 1560 MeV/c. The comparison with the data from 2^ shows 
superscaling for negative values of ip' including ip' < —1, in contrast to the RFC model 
where f{ip') = for ip' < —1. 

(iii) The y-scaling function F{y) is defined in the CDFM on the basis of the RFG relation- 
ships. The calculations of F{y) are performed in the model using three different momentum 
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distributions: from the CDFM, from the ?/-scahng analyses |13l . Il4| and from the theoreti- 
cal approach based on the light-front dynamics method [2^. Comparing the results of the 
calculations for ^He and ^^Fe nuclei with the experimental data, we show the sensitivity of 
the calculated F{y) to the peculiarities of the three n{k)^s in different regions of momenta. 

(iv) An approximate relationship between f{ip') and F{y) is established. It is shown that 
the momentum distribution new for ^^Fe from the y-scaling studies in jisl Q| can describe 
to a large extent the empirical data on f{tp') for q = 1000 MeV/c. We point out that 
the interrelation and the comparison between the results of the and y-scaling analyses 
have to be studied accounting for the correct non-linear dependence of ip' on the ^/-scaling 
variable, which reflects the dependence on the missing energy. 

(v) The regions of momenta in n{k) which are mainly responsible for the description of 
the y- and '?/''-scaling are estimated. It is shown in the present work that the existing data 
on the y- and T/^'-scaling are informative for the momentum distribution n{k) at momenta 
up to ^ 2 ^ 2.5 fm^^. It can be concluded that further experiments are necessary for 
studies of the high-momentum components of the nucleon momentum distribution. 
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